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FIGURE CAPTIONS 



Fig. 1. — The pdf of the density field from the 1.2 Jy iras galaxy catalog (circles connected 
by a solid line) compared with a Gaussian (dotted line) and with the pdf of a model, 
computed from a series of 20 density fields of Gaussian cdm Monte-Carlo simulations 
(triangles). The error bars are the standard deviations over the simulations and measure 
both the effect of the limited volume sampled by iras and the effect of shot noise in the iras 
sampling. Note, however, that the data points are correlated due to smoothing. 

Fig. 2. — The ipdf recovered from the 1.2 Jy iras catalog, compared to the ipdf recovered 
from the Gaussian cdm model, for different values of h. The notation is as in Fig 1. 
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ABSTRACT 

We recover the one-point probability distribution function of the initial 
density fluctuations (ipdf) from the quasi-linear galaxy density field of the 
1.2 Jy IRAS redshift survey smoothed by fO h~^Mpc. The recovery, using the 
laminar, Eulerian, Zel'dovich approximation of Nusser and Dekel, is independent 
of and is relatively insensitive to linear galaxy biasing in the range 0.5 < 6 < 2. 
Errors due to discrete sampling, the limited survey volume, and the method 
of recovering the ipdf are evaluated by comparing the ipdf determined from 
the IRAS data with Monte-Carlo /_Ryi5-like catalogs "observed" from A^-body 
simulations of CDM models with Gaussian initial conditions. Eight sensitive 
statistical tests used in this comparison find the iras ipdf to be consistent with 
Gaussian. We provide observational constraints on possible deviations from 
Gaussianity, which should be obeyed by any theoretical model. 
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1. INTRODUCTION 

The one-point initial probability distribution function (ipdf) is of particular interest 
among the statistics that characterize the primordial density fluctuation field. Whether or 
not it was Gaussian has far-reaching cosmological implications. First, the ipdf provides 
direct information about the origin of the fiuctuations, being independent of the nature of 
the dark matter, unlike the power spectrum, whose evolution during the plasma era depends 
on the physics of the dark matter. Moreover, a known ipdf would place very interesting 
constraints on the value of the cosmological density parameter, 0, based on the observed 
large-scale peculiar velocity field, but independent of galaxy biasing (Nusser & Dekel 1993, 
henceforth ND; Bernardeau et al. 1993). 

The "standard" model, based on infiation, predicts that the initial fiuctuations were 
a Gaussian random field (to all orders of the joint ipdf's, cf., Bardeen et al. 1986), while 
competing models predict various kinds of non-Gaussian fiuctuations. The latter include 
certain infiation models (cf., Kofman et al. 1991), as well as models where the perturbations 
were seeded by cosmic strings (see Bertschinger 1989 for a review), by textures (Turok 
1991) or by non-gravitational explosions (see Ostriker 1988 for a review). 

Even if the ipdf were Gaussian, the present probability distribution function (pdf) is 
expected to be non-Gaussian due to nonlinear effects. The tails of the distribution develop 
positive skewness because high peaks collapse to large densities while the density in voids 
cannot become negative. At moderate amplitudes the fluctuations develop skewness of 
an opposite sign because fluctuations above the mean tend to contract while expanding 
underdense regions occupy more comoving volume as time goes by. 

The detailed quasi-linear evolution of the PDF in an initially Gaussian system has been 
studied using several approximations. For example, according to second-order perturbation 
theory, the third moment of the density distribution grows in proportion to the square of 
the variance (Peebles 1980; Bouchet et al. 1992; Juszkiewicz & Bouchet 1992; Coles & 
Frenk 1992; but see Lahav et al. 1993). Kofman et al. (1993), have found that the density 
PDF actually develops a log-normal shape, as originally proposed by Coles & Jones (1991). 
This is demonstrated in Fig. 1, in which we plot as triangles the mean of the PDF from 
20 A^-body simulations (§4), smoothed by a Gaussian radius of 10 h~^Mpc, and as the 
long-dashed curve a log-normal distribution with the same standard deviation. The error 
bars are the dispersion per bin among the A^-body simulations. 

The present PDF of the density held derived from the 1.2 Jy iras galaxy catalog with 
the same smoothing (§2), is also shown in Fig. 1. The deviations of the present iras pdf 
from Gaussianity agree well with the A^-body simulations. This is a necessary condition for 
Gaussian initial fluctuations, but it is important to realize that the present pdf has only 
limited discriminatory power against non-Gaussianity in the initial conditions. The pdf 
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tends to develop a log-normal shape in a robust way even for certain non-Gaussian initial 
fluctuations {e.g. Weinberg & Cole f992, Fig. f4). 

A more effective strategy to determine the statistical nature of the initial fluctuations 
is to take advantage of the full density field available at the present epoch, trace it back in 
time to recover the initial fiuctuation field (or certain statistics that characterize it), and 
use those to discriminate between theories in the linear regime. One way to accomplish this 
goal is by first recovering the full growing mode of the initial fiuctuation field (e.g., Peebles 
f989, f990; Nusser & Dekel f992; Giavalisco et al. f993), and then computing the ipdf 
from it. Instead, we use below the simple short cut approximation developed by ND in 
which the ipdf is recovered directly from the present quasi-linear density field. We prefer 
not to limit the analysis to certain moments of the distribution, such as the skewness, 
because they might be dominated by the poorly determined tails of the distribution and 
might even diverge. Instead, we recover the full ipdf and only then use several different 
statistics to measure its characteristic features, which may or may not be emphasized by 
the standard moments. 

The shape of the PDF is preserved in the linear regime, so a study of its evolution 
is necessarily nonlinear. On the other hand, we wish to avoid highly nonlinear systems 
in which multi-streaming and mixing obliterate the memory of initial conditions, thus 
rendering impossible the recovery of the ipdf. This restricts us to the quasi-linear regime, 
in which the approximation due to Zel'dovich (1970) is an excellent tool. We follow 
ND in using the Zel'dovich approximation to determine the ipdf from the present day 
velocity field, the latter being calculated from the density field using the phenomenological 
approximation of Nusser et al. (I99I, hereafter NDBB). 

The recovery of the ipdf from a given galaxy-density field requires, in principle, two 
assumptions: (I) a "biasing" relation between galaxy density and mass density and (2) a 
cosmological model, i.e., a value for 0. (The growth of perturbations is insensitive to the 
cosmological constant, e.g., Lahav et al. I99I.) However, in the Zel'dovich approximation 
which we use, the ipdf recovered from a present day density field does not depend on the 
assumed value of (ND, see §3.6 below), unlike the ipdf determined from the present day 
velocity field, which is very sensitive to it (ND). Moreover, the dependence on biasing is 
also quite weak in practice. We are therefore able to recover the ipdf with great confidence, 
with little or no dependence on biasing and on the cosmological model. 

The data are described in §2, the recovery method of the ipdf is detailed in §3, the use 
of A^-body simulations for error analysis is discussed in §4, the recovered ipdf is shown and 
several statistics are computed in §5, and our conclusions are summarized in §6. 
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2. DATA 

The database used here is a complete redshift survey of 5313 galaxies brighter than 
1.2 Jy at 60/im detected by the Infrared Astronomical Satellite [iras)."^ Galaxy candidates 
were chosen from the iras Point Source Catalog, Version 2, (1988) using the selection 
criteria described in Strauss et al. (1990) and Fisher (1992). The data for the brighter half 
of the sample can be found in Strauss et al. (1992). At present, thirty objects (0.5% of the 
sample) remain unobserved. Sky coverage is complete for |6| > 5°, with the exception of a 
small region of the sky which iras failed to survey and regions limited by confusion; our 
sample covers 87.6% of the sky. 

The transformation of a redshift survey to a three-dimensional density map involves a 
number of ingredients. Full details are given by Yahil et al. (1991). The important steps 
are an accurate determination of the selection function, i.e., the probability that a galaxy 
at a given distance be included in the sample, correction for the 12% of the sky not covered 
by the survey, and a self-consistent transformation from the observed redshifts to distances, 
using peculiar velocities deduced from the galaxy distribution itself. This transformation 
from redshift space to real space, in which redshifts are corrected for peculiar velocities, is 
calculated in linear theory, assuming 0°'^/6 = 1. Method 2 of Yahil et al. (1991) is used in 
the present work. 

Our treatment is not completely self-consistent because the transformation from 
redshifts to real distances for the iras galaxies was determined using linear theory (Yahil 
et al. 1991), while the velocity field used for the purpose of the recovery of the ipdf in the 
present paper is calculated from the iras density field using nonlinear theory. However, the 
smoothed iras density field is not very sensitive to this detail. 

The density field of the 1.2 Jy iras galaxies described above was first computed using 
a cloud-in-cell procedure on a cubic grid measuring 240 h~^Mpc on the side, centered on 
the Local Group, with a spacing of 5 h~^Mpc. It was normalized to an average density of 
unity in a sphere of radius 70 h~^Mpc centered on the Local Group. The density was then 
smoothed on the full grid by a Gaussian with radius Rg = 10 h~^Mpc using FFT. In the 
recovery of the ipdf (§3.1-§3.4), however, only a sub-grid of 160 h~^Mpc centered on the 
Local Group was used. In order to reduce errors due to sparse sampling, the final sampling 
of the IPDF (§3.6) was then further restricted to an inner sphere of radius 70 h~^Mpc 
centered on the Local Group. 



■*The Infrared Astronomical Satellite was developed and operated by the U.S. National Aeronautics and 
Space Administration (NASA), the Netherlands Agency for Aerospace Programs (NIVR), and the U.K. 
Science and Engineering Research Council (SERC). 
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RECOVERY METHOD 



The ND recovery method begins with a smoothed density field of galaxies on an 
Eulerian grid at the present epoch and determines the ipdf in the following steps: (a) 
convert the galaxy density field observed at the present epoch to a mass density fiuctuation 
field, 8{x), using an assumed "biasing" relation, (b) solve for the associated peculiar velocity 
field, v{x) using a nonlinear scheme, (c) compute the deformation tensors with respect to 
the Eulerian coordinates at the present epoch at the grid positions, (d) use the Zel'dovich 
approximation to convert them to deformation tensors with respect to the corresponding 
Lagrangian coordinates, (e) obtain the initial linear densities at the above Lagrangian 
coordinates from the Lagrangian deformation tensor, and, (f) use the initial linear densities 
at the Lagrangian coordinates to compute the volume — i.e., mass — weighted distribution of 
density at that epoch. 



If galaxy formation is biased, galaxies need not be faithful tracers of the mass (see 
Dekel & Rees 1987 for a review). There is, however, growing evidence that the galaxy and 
mass densities are strongly correlated (Dekel et al. 1993). Both fields feature the same 
main structures — the Great Attractor, the Perseus-Pisces supercluster and the large void in 
between — and they are consistent with being noisy versions of each other. 

In this paper we assume a deterministic power-law relation between the galaxy and 
mass densities. 



where is the smoothed density fiuctuation field of iras galaxies and 8 is the mass-density 
fiuctuation field. In the limit 8 <^ 1, this model reduces to linear biasing, 8i = h8, which 
is the simplest realization of a linear statistical biasing relation between the variances of 
the fields, as is roughly predicted for linear density peaks in a Gaussian field (Kaiser 1984; 
Bardeen et al. 1986). The power law relation attempts to give a better description of 
biasing, particularly at negative (^'s, where linear biasing can lead to 8 < —1. The nonlinear 
relation between 8j and 8 makes the mean value, (^), deviate slightly from zero (the iras 
density is normalized so [8j) = 0), but this small shift would not affect the ipdf, which is 
computed relative to the mean in the given volume. 

The correlation between the density fields of galaxies and mass has been found by 
Dekel et al. (1993) to be consistent with the above biasing scheme. Motivated by this result, 
we confine the following analysis to 0.5 < & < 2, their 95% confidence limits for < 1.5. 
Our analysis is independent of the value of (§3.6); in practice we use = 1. 



3.1. From Galaxy to Mass density 
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3.2. From Density to Velocity 



In the linear regime the velocity divergence is simply V • w = —af[il)5 where 
/(O) = c/log _D/c/log a 0°-^ is the logarithmic growth rate of the density perturbations 
with respect to the expansion scale, and the gradient is with respect to the comoving 
coordinates. NDBB found empirically that in the quasi-linear regime this divergence is well 
approximated by 



1 + 0.18^ 



Eq. (2) is a Poisson equation in which the density has been replaced by a simple 
function of the density, and since the smoothed velocity field is irrotational for both 
linear and quasi-linear perturbations (e.g, Dekel, Bertschinger, & Faber 1990), it can be 
computed with the usual grid-based FFT techniques. In order to reduce the effect of the 
periodic boundary conditions, we "zero pad" by embedding our 160 h~^Mpc grid in one 
that measures 320 h~^Mpc on the side, with 5 = m the added volume. We do not use the 
IRAS data which we do have in that region, because our smoothing length of 10 h~^Mpc 
becomes smaller there than the mean distance between observed galaxies, and the nonlinear 
correction in Fq. (2) can lead to large errors. 

The accuracy of the velocity field derived using this approximation was tested by NDBB 
and more recently by Mancinelli, Yahil, Ganon & Dekel (1993) in A^-body simulations. 
With Gaussian smoothing of 10 h~^Mpc, the approximation was found to have rms errors 
less than 50 km s~^, compared with the peculiar velocities themselves, which can exceed 
1000 km s-K 



3.3. Eulerian Deformation Tensor 



The Zel'dovich approximation assumes that the comoving Fulerian position s of a 
particle with initial Lagrangian position q can be written as 

x{q,t) = q+D{t)'^{q) , (3) 

where the displacement is separable into a product of the growth function D[t) and a 
time-independent spatial perturbation function, ip. It follows that the velocity is given by 

t> = ax = aD^ . (4) 

We assume further that the fiow is laminar, and that there exists a one-to-one correspondence 
between the Fulerian and Lagrangian positions. The latter can therefore be expressed in 
terms of the Fulerian position and the velocity as 

q{x^t) = X — Dt/){x^t) . (5) 
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The Eulerian deformation tensor, dtpi/dx^^ is readily computed on the grid from the 
velocities determined in §3.2, using Eq. (4) to relate v and t/). The Lagrangian deformation 
tensor, dvi/dqj^ can be obtained from the Eulerian one, using the chain rule, to give 



kj 



D 



dx,, 



dqk 



(6) 



where the expression in the square bracket is the Jacobian matrix obtained from Eq. (5). 



3.4. Initial density 



Eq. (6) can be inverted to yield the full Lagrangian deformation tensor. In fact, the 
initial linear density depends only on its trace 

<^m = -AnTT^ • (7) 

Eor irrotational flow the deformation tensor is symmetric, and its trace therefore equals 
the sum of its eigenvalues. Moreover, Eq. (6) clearly shows that the Eulerian deformation 
tensor and the Jacobian transformation matrix can be simultaneously diagonalized. Hence, 
each eigenvalue of the Lagrangian deformation tensor is given by 

A, = 1^ , (8) 

1 - Doji, 

where [li is the equivalent eigenvalue of the Eulerian deformation tensor at the present 
epoch. The initial density can therefore be written as 

s^^ = -Ih^Y] • (9) 

-Do i 1 ~ -^o/"j 



3.5. From Initial Density to IPDF 



In the last step, we wish to compute the ipdf of 8in, P[5in)^ i.e., the distribution of 
Sin over equal volume elements, which in the linear limit means equal mass elements. But 
the points at which we sample the present density and computed the initial density are 
homogeneously distributed at the present epoch, and not at the initial time. In order to 
correct for this inhomogeneous sampling, we assign each point a weight 1 + S[Xn)^ which, 
by mass conservation, is the relative mass which should be assigned to the Lagrangian 
position q^. In order to avoid errors due to the sparsity of iras sources at large distances, 
we restrict the sample for which we compute the ipdf to those points whose present 
positions, are within a distance of 70 h~^Mpc from the Local Group. 
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Note that there is no guarantee that the mean initial density we compute, ((^m), 
averages to zero. This is because we imposed the condition (5) = (more precisely {5j) = 0) 
on a finite volume at the present epoch, while, in fact, that volume may be a little overdense 
or underdense. We could iteratively correct our normalization until (Si^) = 0. Instead we 
simply shift the Si^s we obtained until their mean equals zero. Dekel (1981, the appendix) 
discusses the justification for this procedure. 

3.6. Q Dependence 

Eq. (9) shows that, apart from overall scale factors, the dependence on can enter the 
initial density Si^ only via the term Dofii. But Eq. (2) depends on only through the scale 
factor a/, and Eq. (4) introduces another factor of 1/aD. Hence, the quantities Dofii are 
given by the product of a function of 5 and Dohf / aDo = 1. The weights used to convert 
the SinS to the final ipdf (§3.5) are functions of ^, so are also independent of 0. The 
dependence on enters only in the transformation of the original iras density from redshift 
to configuration space, and this dependence is quite weak. In our approximation, therefore, 
the IPDF is almost independent of 0. There is a dependency on biasing, introduced in the 
transformation from 5j to ^, but this dependency is also weak (§5). 

4. ERROR ANALYSIS 

The derived ipdf can be in error compared with the universal ipdf for three main 
reasons: (a) the fact that the local volume sampled by the 1.2 Jy iras survey may still 
not be big enough to provide a "fair" sample, (b) shot noise due to the discrete sampling, 
which gets worse at larger distances from the Local Group, and (c) errors in the recovery 
procedure itself. We estimate all these errors simultaneously by computing the ipdf from 
many Monte-Carlo fake catalogs, which were "observed" d la iras from many random 
simulations of a cosmological model. 

The A^-body simulations were run using a particle-mesh code (Bertschinger & Gelb 
1991) with 64"^ grid cells and particles in a periodic cubic box of comoving size 160 h~^Mpc. 
The initial conditions for each simulation were a random Gaussian realization of the 
"standard" cdm spectrum, with = 1 and h = 0.5 (Davis et al. 1985). The simulations 
were halted at a time step when, based on linear growth, the rms density fiuctuation in 
top-hat spheres of radius 8 h~^Mpc equaled a desired value, (Tg,- Since the corresponding 
rms fiuctuation observed for the galaxy distribution is roughly hag, = 1, the value of CTs at 
the final step of the simulation equals the inverse of the biasing parameter, h. We ran 20 
simulations for each of the three values (7^^ = 6 = 0.5, 1,2. 
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We need to construct fake /_Ryi5-like samples from each A^-body simulation, in order 
to compare with the iras reconstruction. We create these samples using the following 
steps: (a) create a smoothed mass density field on a grid by cloud-in-cell, followed by 
Gaussian smoothing with a radius which equals the radius used in the iras reconstruction 
(10 h~^Mpc), (b) bias the mass density field according to Eq. (1), (c) deconvolve the biased 
density field to obtain an unsmoothed galaxy density field, ^ (d) sample the "observed" 
galaxy number density field at each grid point by seeking a random Poisson deviate whose 
expectation value equals the galaxy density obtained in (c), multiplied by ni(/)(r), the 
expected number of galaxies which iras would observe at distance r in a homogeneous 
universe (Yahil et al. 1991), (e) divide the galaxy number density at each grid point by 
ni(f){r) to obtain a new "i^RAS-noisy" distance-independent estimate of the real galaxy 
density. This provides a fake iras catalog, equivalent to the one obtained from the real 
data by the cloud-in-cell procedure. The first smoothing and deconvolution, steps (a)-(c), 
are needed only for biased models, because the biasing has to be applied to the smoothed 
density. On the other hand, the sampling using the iras selection function, steps (d)-(e), 
has to be carried out with the unsmoothed distribution, because the iras simulation so 
obtained will then be smoothed in the ipdf reconstruction, and we have to avoid double 
smoothing. 

In Fig. 1, we show the present PDF of the initially-Gaussian model, with 6=1. The 
PDF is computed for each of the 20 Monte Carlo galaxy density fields. The mean of the 20 
PDF values in each bin is shown as a triangle, and the standard deviation represents the 
error. 

Each of the Monte Carlo galaxy density fields is then used as input for the recovery of 
the IPDF. The mean and standard deviation over the Monte Carlo simulations represent the 
model IPDF and the uncertainties about it. Eig. 2 shows the ipdf recovered from the model, 
and the associated errors, for the three different values of the biasing factor, h = 0.5, 1,2, 
which are the "true" values used in the simulations themselves. 

5. RESULTS 

The IRAS galaxy density field is fed into the recovery procedure, which yields the ipdf 
shown as solid curves in Eig. 2 for the three assumed values of h. We can see that the ipdf 
is quite insensitive to the assumed value of h in the relevant range [0.5, 2]. This insensitivity 
to the assumed 6, allows us, without loss of generality, to compare the iras ipdf only to 
model iPDF^s in which the recovery assumed the "true" value of b. We clearly see in Eig. 2 
that the iras ipdf is consistent with a Gaussian, for all the three values of b tested. 



®The deconvolution is unstable if the broadening function in fc-space is allowed to go to zero at high k. 
We modify the fc-space Gaussian by multiplying it by 0.99 and adding a 0.01 white background. 
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We next wish to quantify the agreement, or discrepancy, between the data and the 
Gaussian modeL Quantitative hmits on the possible deviation from Gaussianity, by several 
different statistics, would be particularly useful for testing non-Gaussian theoretical models. 
The counts in the different bins that construct the ipdf are correlated in a complex 
fashion because the fields were smoothed over a large scale at the present epoch before the 
reconstruction was begun. In practice, such complex correlations can only be analyzed by 
using statistics that measure certain properties of the ipdf, and evaluate the distributions 
of these statistics using Monte Carlo simulations. We determine the standard deviation of 
each statistic over the Monte-Carlo catalogs, and assume a Gaussian distribution for this 
statistic when we express the confidence level of rejection in terms of numbers of sigma. 

The results of several statistics are summarized in Table 1, which lists for each value of 
h and for each statistic the mean value and the standard deviation in the Gaussian model, 
the value measured from the iras data, and its deviation from the value obtained for the 
Gaussian model in units of the standard deviation. 



First listed are [i and cr, the mean and standard deviation of Si^. Then comes W , a 
standard statistic which has been optimized to measure deviations from Gaussianity when 
the distribution is sampled by a small number of independent points. Given the initial 
densities at grid points, z = l,n, the statistic is defined by 



where [n/2] indicates the integer part of n/2, and the coefficients are chosen to maximize 
the discriminatory power against Gaussianity (Shapiro & Wilk 1965). The coefficients are 
computed according to the prescription described by these authors. The initial density at 
each grid point was weighted by the present density there. 

There follows a x^-type statistic, which measures a vertical deviation between the data 
IPDF and the model ipdf: 



where i runs through all the bins and (Ti is the total error of Pi in that bin. This is, of 



The rest are higher moments of the distribution. Let x = (^j„ — fj,)/crs^^. By definition 
(x) = and (x^) = 1. The skewness and kurtosis are commonly defined by 



[Table 1] 




(10) 





course, not a proper x 



^ because the data points are correlated. 





such that they both vanish for a Gaussian distribution. We see from Table 1 that these 
moments carry large random errors, on the order of a few tenths. This is because they 
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are dominated by the poorly determined tails of the distribution. For measuring similar 
properties of the ipdf but with less weight assigned to the tails, ND suggested replacing 
the skewness and kurtosis by 



The former, ranging from —1 to +1, measures mirror asymmetry about zero. The latter, 
ranging from — ^2/7r to 1 — ^2/7r, measures symmetric deviations from a Gaussian; it is 
negative when the PDF is dominated by excessive tails or by a shallow center. We see in 
Table 1 that the errors of these moments are much better behaved; their random errors are 
only on the order of a few percents. 

As expected on the basis of a visual inspection of Fig. 2, the Gaussian model cannot 
be ruled out by the data using any of the statistics. The most discriminatory among the 
statistics seems to be W, by which the model with 6 = 1 is "rejected" at the 1.85-sigma 
level. In all other cases the iras datum lies well within the very likely region of the 
distribution of the statistic. The ipdf is thus fully consistent with a normal curve. This 
conclusion is insensitive to the value of 6, at least in the range [0.5,2]. 

In order to quote n-sigma limits on the allowed deviation from Gaussianity based on 
any statistic s, one can compute from Table 1 the quantity (s — (s)) it ncr^, where (s) and 
(Tg are the model Monte-Carlo mean and standard deviation of s. For example. For a 6 = 1 



model we find the 3-sigma constraints, 0.87 <W< 1.03, < 1-2, -0.18 < Sr < 0.11, 
-0.205 < Kr < -0.097, -0.65 < 5" < 0.36, and -0.82 < K < 0.62. 



Based on the 1.2 Jy iras redshift catalog, and for biasing-parameter values in the 
range of interest, 0.5 < & < 2, the density ipdf is fully consistent with a Gaussian. Any 
possible deviations from a Gaussian are confined by several sensitive statistics to a small 
permitted range; these statistics could also be useful when trying to evaluate the likelihood 
of non-Gaussian models. The method used here, based on a recovery of the ipdf, is more 
sensitive than tests using the present PDF. 

One limitation of the current results is the model dependence of the error analysis, in 
the Monte-Carlo analysis, we had to use a Gaussian model with a specific power spectrum. 
Although the power spectrum is expected to have a negligible effect on the ipdf itself, it 
should affect, for example, the estimated errors due to the limited volume sampled. The 
CDM spectrum used here is an appropriate choice because it is consistent with the observed 
power spectrum over the scales used in this investigation (Vogeley et al. 1992; Fisher et al. 
1993; Baugh & Ffstathiou 1993; Kolatt & Dekel 1994). It would nevertheless be worthwhile 
to test the method and estimate the errors using other models with different spectra. 




(13) 



6. 



DISCUSSION AND CONCLUSIONS 
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Another limitation is that the constraints on possible deviations from a Gaussian ipdf 
were derived based on Monte-Carlo simulations of a Gaussian model. The appropriate 
procedure for testing a specific non-Gaussian model is to estimate the distribution of each 
statistic by Monte-Carlo simulations perturbed from that specific model. The constraints 
provided by Table 1 are therefore reliable only for models that do not deviate drastically 
from Gaussian. Any testing of a strongly non-Gaussian model would require a repetition of 
the full Monte-Carlo analysis using the ND recovery technique. 

The Gaussian ipdf found here from iras adds a crucial ingredient to the framework 
discussed in ND. Quasi-Linear gravity is used to relate two observables: (1) a smoothed 
velocity field and (2) a smoothed galaxy-density field, to three fundamental theoretical 
ingredients: (a) the cosmological model via 0, (b) the relation between galaxy and mass 
density via the biasing parameter 6, and (c) the initial fiuctuations via their ipdf. Given 
a galaxy-density field, one can recover the ipdf independently of and with only weak 
dependence on h. Given a velocity field, one can assume a value for and obtain the ipdf, 
independently of biasing. Given both fields, one can simultaneously constrain both and 
the IPDF. 

The contribution of the current analysis is to establish the Gaussian nature of the 
initial fiuctuations, essentially independently of the assumed values of and h. This, 
in turn, argues that is about unity or larger, because ND showed that the potfnt 
velocity field is consistent with a Gaussian ipdf if and only if has such a large value. 
For example, values of < 0.3 were ruled out at the 4-6 sigma level by several different 
statistics. Similar, though somewhat weaker conclusions, were reached by Bernardeau et al. 
(1993) based on the skewness of V • w in the potfnt reconstruction and the assumption of 
Gaussian initial fluctuations. 

Another two of the theoretical ingredients, and 6, can be determined from a direct 
quasi-linear comparison of the POTSJVT-mass and i^RAS-galaxy density flelds. The strongest 
result of Dekel et al. (1993) is 0°'^/6 = 1.28^0^59 at the 95% confldence limit. Separate 
nonlinear constraints on and 6, which are somewhat weaker, are > 0.46 for h > 0.5 
and b = 0.7t.o'^ for = 1, again at the 95% confldence level. These constraints, as well as 
the ones from the potfnt reconstruction of the ipdf, are expected to tighten shortly with 
forthcoming new velocity data. 

Once the analysis is sufflciently nonlinear, we expect only one simultaneous solution for 
all three theoretical ingredients. The solution we have found so far from the above-mentioned 
investigations combined, i.e., ^ 1, 6/ ^ 1, and Gaussian ipdf, are all consistent. The 
constraints therefore seem to be closing in on a "standard" model of an Einstein-de Sitter 
flat universe with Gaussian initial fluctuations and iras galaxies roughly tracing mass 
on scales ~ 10 h~^Mpc. Low values of 0, with or without a cosmological constant, or 
non-Gaussian initial fluctuations, cannot be easily invoked when trying to explain the 
remaining puzzles of the formation of large-scale structure in the universe. 
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